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Abstract 

O . 

{Nj , This paper is concerned with the quantum-mechanical decay of a Schwarzschild-like 

black hole, formed by gravitational collapse, into almost-flat space-time and weak radi- 
; ation at a very late time. We evaluate quantum amplitudes (not just probabilities) for 
transitions from initial to final states. This quantum description shows that no informa- 



(N 



tion is lost in collapse to a black hole. Boundary data for the gravitational field and (in 
this paper) a scalar field are posed on an initial space-like hypersurface E/ and a final 
surface Ej? . These asymptotically-flat 3-surfaces are separated by a Lorentzian proper- 
time interval T (typically very large) , as measured at spatial infinity. The boundary-value 
problem is made well-posed, both classically and quantum-mechanically, by a rotation of 
T into the lower-half complex plane: T — > \T\ exp(-iO), with < 9 < n/2 . This corre- 
sponds to Feynman's + ie prescription. We consider the classical boundary-value problem 

■ (2) 

and calculate the second-variation classical Lorentzian action S„J aa as a functional of the 
| boundary data. Following Feynman, the Lorentzian quantum amplitude is recovered in 
the limit 9 — > 0+ from the well-defined complex-T amplitude. Dirac's canonical ap- 
proach to the quantisation of constrained systems shows that, for locally-super symmetric 
theories of gravity, the amplitude is exactly semi-classical, namely, exp(£S^j2g S ) for weak 
perturbations, apart from delta-functionals of the supersymmetry constraints. We treat 
such quantum amplitudes for weak scalar-field configurations on Sp, taking (for simplic- 
ity) the weak final gravitational field to be spherically symmetric. The treatment involves 
^ ■ adiabatic solutions of the scalar wave equation. This considerably extends work reported 
in previous papers, by giving explicit expressions for the real and imaginary parts of such 
quantum amplitudes. 

PACS numbers: 04.70.Dy, 04.60.-m 

1. Introduction 

This paper is concerned with the simplest example of quantum radiation following 
gravitational collapse to a black hole, in which the Lagrangian contains only Einstein 
gravity and a minimally-coupled massless scalar field (p . As with emitted particles of spin 
s > [1-3] , so here we find that the emission process with s = may be described in terms 
of quantum amplitudes (not just probabilities) [1,4-6]. An effective-field approach to this 
question has also been described — see [7,8]. 

We compute such amplitudes; the basic quantity is the amplitude to go from data on 
an initial asymptotically-flat space-like hypersurface E/ to data on a final surface Ei? . 
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The proper-time separation between E j and E_p , as measured at spatial infinity, is T . The 
4-dimensional space-time metric is denoted by (//, v = 0,1,2,3). The gravitational 
data to be specified on E/ and on E^ (regarded as surfaces of constant coordinate x ) are 
the intrinsic positive-definite spatial 3-metric hij = gij on each 3-surface, written as hiji 
and hijF ■ The scalar data are taken to be the scalar field 0/ on E/ and 4>f on Hp ■ 

In field theory, quantum calculations are usually carried out about some classical 
solution. The above boundary data with T real, however, are clearly ill-suited to such a 
classical boundary-value calculation, since classical hyperbolic (wave-like) equations are, in 
general, badly posed when subject to data on an outer boundary [9] . (The natural arena for 
hyperbolic equations is in the Cauchy evolution problem [10,11].) A simple example of such 
badly-posed behaviour is given in [5] (see also Sec. 2. 4 of [12]), which treats the 1-space, 
1-time weak-field analogue of the boundary-value problem of the previous paragraph. In 
this example, there is only a real massless scalar field in flat Minkowski space-time. The 
roles of Ej and E^ are taken by the pair of lines {t = , T ; — oo < x < oo}, and, 
for simplicity, data for are chosen with <p{t = 0,x) = 0, <p(t = T, x) = (f>i(x), where 
4>i is a specified function of rapid decrease as \x\ — > oo . One can see from [5] that 
this boundary-value problem becomes well-posed, once T is rotated into the lower-half 
complex plane: 

T — ► |T|exp(-z0) (0 < 9 < tt/2). (1.1) 

However, the limit 9 — > is somewhat singular. This fits in with the behaviour that one 
might expect: for 9 exactly equal to tt/2 , one has a real elliptic boundary-value problem 
for the Laplace equation, which is known to be well-posed. In between, for < 9 < n/2 , 
the differential equation has complex coefficients but the boundary-value problem is still 
well-posed, with existence, uniqueness and analyticity of the solutions. We are, in fact, 
viewing a strongly elliptic boundary-value problem when < 9 < tv/2 [13]. For 9 = , one 
is, as above, trying to fit a wave-like system into a boundary-value problem. 

Correspondingly, in our field-theoretic black-hole evaporation problem, we make the 
same complex rotation (1.1) of T . Weak-field analyses suggest that, with such a complex 
rotation, the classical Dirichlet boundary-value problem for a perturbative scalar field 
4> (and also for a linearised gravitational field Sg^) become well-posed. Hence, for an 
evolution of the black hole which does not deviate too far from the spherical, one expects 
to be able to study a semi-classical expansion of the amplitude, provided Im(T) < . 
Once one has found the quantum amplitude for a given value of 9 , with < 9 < n/2 , one 
can rotate 9 back towards zero, following Feynman's + ie prescription [14,15]. In the 
context of scalar-field evaporation from black holes, this procedure will be treated in Sec. 8 
below. The slight complexification of T induces an imaginary part in the total Lorentzian 
action, which is crucial in computing the appropriate quantum amplitude. Conversely, even 
for fairly small values of 9 , solution of the rotated classical boundary-value problem is 
expected to smooth any variations or oscillations of the boundary data, as one moves from 
the boundary E/ or E^ into the interior by a few multiples of the relevant wavelength. 

To fix one's physical intuition, imagine that the initial distribution of scalar field 0/ 
on E/ is nearly spherically-symmetric and extremely diffuse, with almost all of the mass 
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distributed over radii much greater than the 'Schwarzschild radius' 2M . Here, M is the 
ADM (Arnowitt-Deser-Misner) mass, which is defined in terms of the rate of fall-off of 
the initial data, at large radii on a given spacelike hypersurface [16]. The initial 3-metric 
hij on Ej will be almost spherically symmetric, and will vary very slowly with radius. The 
final surface E p will be chosen at very late T , so as to register all the evaporated radiation. 
The total ADM mass on the final surface must also equal Mq , since otherwise the classical 
boundary-value problem with a given (finite) time-interval T at spatial infinity will have 
no solution [17,18]. As regards the classical solution in the interior (0 < t < T), the 
geometry is well approximated at late times by the radiating Schwarzschild-like Vaidya 
metric [19-21]. The classical scalar-field solution will depend on the enormous amount of 
detail which, in general, is present in the prescribed final data ^ on Sp . 

In Sec. 2, we comment on the need to work, not just with the Einstein/scalar La- 
grangian, but with supergravity or with gauge-invariant supergravity models which in- 
clude supermatter, in order that quantum amplitudes should be meaningful. Indeed, in 
the locally-supersymmetric case, for a large class of models, the amplitude above turns 
out to be exactly semi-classical, in a certain sense — see Eq.(2.2) below [12,22-24]. In 
Sec. 3, we discuss the 'background' nearly-spherically-symmetric 4-metric 7^ , needed in 
a self-consistent treatment of the classical field equations. In particular, the Einstein field 
equations give, at lowest order, a 'source' for 7 MI/ which includes the energy-momentum 
tensor of the scalar field <f> , together with a source quadratic in graviton perturbations 
(as well as corresponding sources for any other matter fields present). Sec. 4 treats the 
decomposition of scalar perturbations in spherical harmonics, assuming that the back- 
ground 4-metric 7^ is spherically symmetric. In Sec. 5, we describe the classical action 
functionals S c \ a ss (Lorentzian action) or I c i a ss (Euclidean or Riemannian action), related 
by iS c \ ass = — /class 7 for the Einstein/scalar system. Because S c \ ass or I c i a ss are evaluated 
at a solution of the classical field equations, they reduce to a sum of boundary terms. Sec. 6 
treats the adiabatic radial functions for evolution of the (linearised) scalar field, based on 
the general treatment in Sec. 4. For adiabatic perturbations (high frequencies), the time- 
dependence is approximately harmonic and can be factored out, leading to a second-order 
radial equation for given frequency u and angular quantum numbers I , m . In Sec. 6, we 
also describe the 'coordinates' which are most convenient in specifying the final data ^on 
Ef . A suitable basis of radial eigenfunctions on the final surface E^ is discussed in Sec. 7. 
The analytic continuation process, in which a Lorentzian quantum amplitude is derived 
from the limit 9 — > 0+ of the amplitude for complex time-separation T, as in Eq.(l.l) 

[14,15], is treated in Sec. 8. By these methods, one can, if desired, evaluate both the real 

(2) 

and imaginary parts of the lowest-order perturbative classical action S^iass , and hence of 

( 2) 

the semi-classical amplitude exp(i5' i : lass ). Sec. 9 contains the Conclusion. 

2. The quantum amplitude for bosonic boundary data 

Consider, at present, the 'Euclidean' quantum amplitude to go between prescribed ini- 
tial and final purely bosonic (gravitational and matter) data, on a pair of 3-surfaces, each 
'topologically' (diffeomorphically) M 3 and each carrying an asymptotically-flat 3-metric. 
One further needs to specify the proper (Euclidean) distance r , measured orthogonally 
between the two surfaces at spatial infinity. From one (purely formal) point of view, this 
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amplitude can be regarded as given by a Feynman path integral over all Riemannian infill- 
ing 4-geometries, together with any other fields, each such configuration being weighted 
by exp(— I), where / is the 'Euclidean action' of the configuration. If this definition were 
meaningful, one would expect that the resulting 'Euclidean' quantum amplitude would 
have the semi-classical form 

Amplitude ~ (A + HA 1 + h 2 A 2 + . . . ) exp(- I B /h), (2.1) 

giving an asymptotic expansion in the limit that (Ib / K) — > . Here, Ib is the classical 
'Euclidean action' of a Riemannian solution of the coupled Einstein and matter classical 
field equations, subject to the boundary conditions. For simplicity, we assume that there 
is a unique classical solution, up to gauge and coordinate transformations. It is quite 
feasible, though, in certain theories and for certain boundary data, to have instead (say) 
a complex-conjugate pair of classical solutions [25]. 

The classical action Ib and loop terms Aq , A\ , A 2 , ... depend in principle on the 
boundary data. In the case of matter coupled to Einstein gravity, each of Ib , Aq , A\ , 
will also obey differential constraints connected with the local coordinate invariance of the 
theory and with any other local invariances such as gauge invariance (if appropriate) [12,26]. 
This follows from the Dirac canonical approach to the quantisation of theories with local 
or gauge-like invariances [12,26]; the 'differential' Dirac approach is dual to the 'integral' 
Feynman approach. The most striking results of the Dirac approach arise when the invari- 
ance properties include local supersymmetry [12,22,23]. For the present boundary-value 
problem, the Dirac approach is much more powerful than the path-integral approach. Fol- 
lowing the Dirac approach, and given local supersymmetry, the semi-classical expansion 
(2.1) simplifies. For example, in N = 1 supergravity, one has, from the Dirac constrained- 
quantisation approach [12,24]: 

Amplitude = A exp(- Is/ft), (2.2) 

apart from factors which are delta-functionals of the classical supersymmetry constraints 
[12,24] on the bounding surfaces. In this theory, the one-loop factor Aq is in fact a 
constant. In Eq.(2.2), Ib denotes the classical action, including both bosonic and fermionic 
contributions. Related semi-classical behaviour holds for A = 1 supergravity coupled to 
gauge-invariant supermatter [22,24,27]. 

In the case (2.2) of N = 1 supergravity, the classical action is all that is needed for the 
quantum calculation. A corresponding situation arises with ultra-high-energy collisions, 
whether between black holes [28] , in particle scattering [29] , or in string theory [30] . 

In the asymptotically-flat, spatially-M 3 context appropriate here, the purely Rieman- 
nian case corresponds, in Lorentzian-time language, to a time-separation at spatial infinity 
of the rotated form T = — ir , where r is the (positive) imaginary-time separation defined 
above. If the four-dimensional classical bosonic part of the solution is to be real, then 
the bosonic boundary data must be chosen real. Following the standard route, one should 
study the (now complex) amplitude (2.1) or (2.2), as a function of the angle 9 of Eq.(l.l), 
where 9 is rotated from 9 = n/2 to 9 = + e (e > 0), with 

T = rexp(-z0). (2.3) 
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Provided that there is a (complex) classical solution to the Dirichlet problem, which varies 
smoothly with 6 , (e < 9 < 7r/2), the expression (2.1) or (2.2) should continue to give the 
quantum amplitude. In particular, this would occur if strong ellipticity [13] held for the 
coupled Einstein/bosonic-matter field equations, up to gauge. 

In the present paper, we assume that the Lagrangian is indeed invariant under local 
supersymmetry, as, for example, in gauge-invariant N = 1 supergravity [27], but (for 
simplicity) that only the gravitational and scalar data are present in the boundary data 
and for the classical solution. 

3. The approximate 4— dimensional metric 

As above, the classical background bosonic fields will here, for simplicity, be taken to 
be simply the metric g^ v and massless scalar field (p . In other work, we study cases in which 
fields with different spins s = \ , 1 or | are included as perturbations of a spherically- 
symmetric background solution [2,3,31]. The classical solutions (g^ , 4>) of the coupled 
Einstein/scalar field equations below are taken to have a 'background' time-dependent 
spherically-symmetric part (7^ , $), together with a 'small' perturbative part (h^ v , </> pe rt)- 
The perturbative fields and <^ pe rt , which live on the spherically-symmetric background 
4-geometry with metric 7^ , can, as usual, be expanded out in terms of sums over tensor 
(spin-2), vector (spin-1) and scalar harmonics [32,33]. Each harmonic is weighted by a 
function of the Riemannian time and radial coordinates (r, r). 

The Einstein field equations read 

Gp,v = R^u — -^Rg^,u = 8tvT^ u , (3.1) 

where R^ u denotes the Ricci tensor, R the Ricci scalar and T^ v the energy-momentum 
tensor. For a real massless scalar field </> , one has 

Tpv = JAt 0,„ - ^g^u (0,a0,/3 9 al3 )- (3-2) 

The gravitational field equations further imply the scalar field equation (the Laplace- 
Beltrami equation [34]): 

g^u) = 0, (3.3) 

where g denotes det(g^), and (for the moment) we assume that the 4-metric g^ v is real 
Riemannian, whence g > . 

The corresponding variational principle involves an action functional of the form [12] 
I = — — J d 4 x g? R + — J d 4 x g? (V</>) 2 + boundary contributions. (3.4) 

The appropriate boundary terms will be discussed in Sec. 5. 

In the Riemannian case [35], the 'large' or 'background' 4-metric can be put in the 
form: 

ds 2 = e b dr 2 + e a dr 2 + r 2 (d9 2 + sin 2 6> # 2 ) , (3.5) 
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where 



b = b(r,r) ; a = a(r, r). 



(3.6) 



If the gravitational field were exactly spherically symmetric, as in Eq.(3.5), and if the scalar 
field were also spherically symmetric, of the form 0(r, r), then the Riemannian spherically- 
symmetric scalar and Einstein field equations would hold [35]. The scalar field equation 
reads: 

+ e b ~ a <(>" + I (a - b) + r" 1 e b ~ a (l + e a )0' = , (3.7) 

where ( ' ) denotes d( )/dr and ( )' denotes d( )/dr . Together with Eq.(3.7), a slightly 
redundant set of gravitational field equations is given by: 

a' = - 47ir (e a - b 4> 2 - </>' 2 ) + r' 1 (l -e a ), (3.8) 
b' = - Anr (e a ~ h 4> 2 - </>' 2 ) - r' 1 (l -e a ), (3.9) 
a = 8yrr00', (3.10) 

a + e b ~ a b" + \ (a - b) a - r~ l e b ~ a (l - e a ) (b' + 2r~ l ) = 8n (0 2 + e b ~ a <j)' 2 ). (3.11) 

The metric and classical field equations in Lorentzian signature [36], or for certain 
types of complex metrics, can be derived from the above by the formal replacement 

* = re"™, (3.12) 

where a is independent of 4-dimensional position, and should be rotated from to n/2 . 

In the bosonic black-hole evaporation problem, the classical Riemannian metric and 
scalar field will not be exactly spherically symmetric; similarly for any non-zero spin-| and 
spin-| classical (odd Grassmann-algebra-valued [12]) fermionic solutions in any locally- 
supersymmetric generalisation [27]. In particle language, rather than the field language 
mostly used in this paper, huge numbers of gravitons and scalar particles will continually 
be given off by the black hole (together with any fermions allowed by the model), leading 
effectively to a stochastic distribution, in which, for any given spin s , the field fluctuates 
around a spherically-symmetric reference field. 

Consider, for example, gravitational and scalar perturbations about a Riemannian 
spherically-symmetric reference 4-metric 7 MI/ and reference scalar field $ . In perturba- 
tion theory in general relativity [37] , one considers a one-parameter (or many-parameter) 
family of 4-metrics, here given by the asymptotic expansion 

9^(x, e) ~ 7/iI/ (x) + e h^x) + e 2 h^x) + . . . , (3.13) 

where h^ 1 ^ is the first-order metric perturbation, h^ 2 \ v the second-order perturbation, etc. 
Throughout, the superscript (°) will refer to the background, while ^ denotes terms linear 
in perturbations, etc. Indices are to be raised and lowered using the background metric 
Covariant derivatives with respect to the background geometry are denoted 
either by a semi-colon ; a or equivalent ly by V a • 
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Analogously, we split a real massless scalar field <j) into a spherically-symmetric back- 
ground piece $(r, r) and a (non-spherical) perturbation: 

(j>(x,e) ~ $(r,r) + e</> (1) (z) + e 2 (2) (a;) + ... . (3.14) 

The spherically-symmetric background part $ will be non-zero if the background scalar 
data at early and late Euclidean times r contain a non-trivial spherically-symmetric 
component. The perturbation fields (j)^(x), . . . will, in general, contain all non- 

spherical angular harmonics. These fields must be chosen such that the entire coupled 
Einstein/scalar system satisfies the classical field equations, as well as agreeing with the 
prescribed small non-spherical perturbations in the initial and final data, both gravita- 
tional and scalar. The effective energy-momentum source for the spherically-symmetric 
part 7^ of the metric includes contributions formed quadratically from the non-spherical 

gravitational and scalar expressions h { % and <j>M — see Eqs. (3.27-29) below for further 
detail. 

In the simplest case, one can restrict attention to the exactly spherically-symmetric 
Riemannian model of Eqs. (3. 5-11). The background metric 7 Mi , corresponds to the metric 
(3.5,6), with respect to suitable coordinates, and <3>(t, r) above corresponds to 0(r,r) of 
those equations. This Riemannian boundary-value problem, involving a system of coupled 
partial differential equations in two variables (r, r), has been studied numerically in [35] for 
particular choices of boundary data, and is currently being investigated in greater detail 
[38]. 

In contrast, the Lorentzian-signature version of the spherically-symmetric classical 
Einstein/scalar system must be studied as an initial-value evolution problem, in order to 
be well posed [36]. One conceivable initial profile for the scalar field, which has been much 
studied in the Lorentzian-signature numerical problem [39,40], is an ingoing 'Gaussian' 
shell of scalar radiation. To define such initial data, work in a nearly-flat space-time at 
very early times (large negative Lorentzian time-coordinate t). Define an advanced null 
coordinate 

v = t + r. (3.15) 
The incoming 'Gaussian' shell is asymptotically, at early times, of the form: 



$(t,r) 




(3.16) 



where fo , k , d , r$ and A are all positive real parameters. The radial extent, Lq , of the 
'Gaussian' is given by Lq ~ A . The numerical evolution of such initial data provides 
a model of spherical collapse. In particular, two main qualitatively different regimes of 
initial data can be distinguished. First, if Lq or A is too large, then the initial data 
are 'diffuse', there is little self-interaction, and the incoming scalar profiles pass more or 
less straight through each other, leaving behind nearly-flat space-time plus small per- 
turbations. Second, if A (or Lq) is less than a certain critical value, the interaction is 
sufficiently non-linear that a black hole forms. 
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Returning to the Riemannian or the complex case, one can expand out the Einstein 
field equations (3.1,2) in powers of e . At lowest order (e°), one has the background Einstein 
and scalar-field equations 

R ( °l- \R {0) 1„ = SttT^ , (3.17) 
7^$;^ = 0, (3.18) 

where i?^„ denotes the Ricci tensor and the Ricci scalar of the background geometry 
7^ . Further, 

T% = ^ $,„ - 1 lixv ($, Q d> )/3 7 «/3) (3.19) 

denotes the background spherically-symmetric energy-momentum tensor. These field 
equations are equivalent to Eqs.(3.7-ll), when coordinates are taken as in Eqs.(3.5,6). 

The linearised (e 1 ) part of the Einstein equations reads (see Section 35.13 of [16]): 



h {1 l^ - 2 >* ( W> - 2 * (0 U h{1) ™ - 2 R{0 \ hil l) * 



(3.20) 



Here, is defined by 

h (1 l = - l^h^, (3.21) 

where 

= feWf . (3.22) 

Also, R^a^isa denotes the Riemann tensor of the background geometry 7aii/ . Further, 
denotes the linearisation or 0(e 1 ) part of the energy-momentum tensor T^ v (x, e). 
Explicitly, 

T { %„ = 2V ( ^ 1 )V, ) $ - 7m-V q $ V>W 

1 (3 23) 

+ - (7m^ ^ W " Y P ) V p $ . 

The linearised Einstein equations (3.20-23) are most easily studied in a 'linearised har- 
monic gauge' [16,41] in which, by an infinitesimal coordinate transformation, one arranges 
that 

h (1 U a = (3.24) 

everywhere. Since the gravitational background 7mi/ is spherically symmetric, the lin- 
earised Einstein equations (3.20-23) can be further decomposed into three independent 
sets of equations. These describe, repectively, scalar (spin-0) perturbations associated 
with matter-density changes (T ^ r ), vector (spin-1) perturbations associated with matter- 
velocity changes (T^J, and gravitational radiation (spin-2) associated with anisotropic 
stresses (T^-) [42]. These equations and their solutions are described further in [2]. 



The linearised (e 1 ) part of the scalar field equation (3.3) yields 

7^0 (1 ;^ " (h^$,„) = 0. (3.25) 

The linearised Einstein and linearised scalar-field equations (3.20-23,25) are coupled. 

At order e 2 , the gravitational field equations give (implicitly, below) the second-order 

(2) (2) 

contribution G ^ to the Einstein tensor G^ u . Note that G ^ includes a well-known 

contribution quadratic in the first-order perturbations h \ v and their derivatives — see 
Eq. (35.58b) of [16]. This part represents an effective energy-momentum-stress density due 

(2) 

to the gravitational perturbations, including gravitons. G ^ also contains contributions at 
quadratic order, formed from the background $ and the linearised <f>^ or their derivatives, 
together with 7 Mi , and h^\ v . These parts represent the 0(e 2 ) contribution of the scalar- 
field energy-momentum tensor T^ v of Eq.(3.2). 

Explicitly, one finds, after a lengthy calculation [1], that the Einstein equations, to 
quadratic order in perturbations, read 

G% = 8n r(%, + 8tt T% + 8tt - G% . (3.25) 

T { % = V/v/)-ivfV/)v/ 

+ (7m^ - h% r P ) V p 0« (3.27) 

+ \ {h {1 U - 1[xv h^ a h^ p ) V CT $ V p $ 

+ -h^\ R<-°\ h^--h^ h {1 \ i*<°> - Srr h {1 \ p 
- 4tt 7 ^ (2 V CT </> (1) V p $ + </> (1) V CT V CT (1) - /i (1)fjp " V p $ V^) 

(3.28) 

where the explicit form of will not be used here. Also, G ^ is defined implicitly by 
Eq.(3.20). The above expressions are needed particularly in studying the Vaidya metric 
[19], which, as shown in [21], describes approximately the late-time region of the geometry 
following gravitational collapse to a black hole, containing a nearly-steady outgoing flux of 
radiation. The Einstein field equations, averaged over small regions, give the contribution 
of massless scalar particles, gravitons, etc., to the nearly-isotropic flux. 

Physically, for the Riemannian or complex boundary-value problem discussed in Secs.l 
and 2, the 0(e) perturbations in the 4-metric g^ u and scalar field 0, relative to the 
spherically-symmetric background solution (7^ , <&), should arise classically from 0(e) 



Here, 



and 
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perturbations away from spherical symmetry in the boundary data , <f> (or d&/dn) at 
the initial and final surface. Moreover, as mentioned earlier, provided that the perturbed 
boundary data contain numerous high harmonics, the 4-dimensional perturbations in the 
interior would be expected effectively to have a stochastic nature. When averaged over a 
number of wavelengths, the effective perturbative energy-momentum tensor above, T EF 1 F , 
will yield a spherically-symmetric smoothed-out quantity < T EFF V > [43,44]. In a locally- 
super symmetric version of this theory, the energy-momentum tensor due to the spin-| 
and spin-| fields will also contribute to < T EF I F >. In particular, this averaged form 
of T EFF will account for the gradual loss of mass by radiation of a black hole in the 
nearly-Lorentzian sector (that is, in the case of a time-interval at infinity of the form 
T = t exp(— i6), where r is large and positive, 6 = e is small and positive). Although < 
T EF l F > is small (being of order e 2 ), its effects on the black-hole geometry, including those 
on the mass, will build up in a secular fashion, over a time-scale of order 0(e -2 ). Such 
secular behaviour appears often in perturbation problems [45,46] - for example, in the 
familiar treatment of the perihelion precession of nearly-circular orbits in the Schwarzschild 
geometry [47] . In our boundary-value problem, whether regarded as classical or quantum, 
the initial boundary data will be spread over a 'background' extent of 0(1) in the radial 
coordinate r on the initial surface Ej . But, corresponding to the 0(e~ 2 ) time-scale for 
the black hole to radiate, the final data on E^ will be spread over a radial-coordinate scale 
of 0(e 2 ). Thus, even the classical boundary-value problem here is an example of singular 
perturbation theory [45,46]. 

The standard treatment of high-frequency averaging in general relativity was given 
by Brill and Hartle [43] and by Isaacson [44]. Let < > denote an average over a time 
To much longer than typical wave periods, together with a spatial average over several 
wavelengths A . Then: 

< 9^ > = Ivlu , <</>>=$, (3.29) 

< </> (1) > = , < h% v > = , (3.30) 

< d a hS% > = , < d a d p h% u > = . (3.31) 

Indeed, 

< CW > = c(o) 5 (3.32) 

for any background quantity . Rules for manipulating these averages in the high- 
frequency aproximation are set out in [44]. Under integrals, the average of total divergences 
can be neglected. For example, 

< /* (1)Q / h(1 l;a > = - < h {1 l > . (3.33) 

Further, covariant derivatives commute for high-frequency waves. The rules (3.29-33) 
imply that 

<T%> = 0, (3.34) 

< T^l > = < V^ 1 ) - \ llw V a </>« > . (3.35) 
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We can now rewrite the background field equations (3.7-11) in a form smoothed out 
by averaging over a number of wavelengths of the scalar and gravitational perturbations 
[43,44]. The equation which includes the quadratic-order contribution of the perturbations 
as a source for the background geometry reads: 

G^Ul) = Stt T<%, + Stt e 2 (< T<%, > + < >). (3.37) 

- (2) 

The terms in this equation vary over length-scales ^> A . The 'source equation' for h ^ , 
analogous to Eq.(3.20) for , is 

G {1 l(l,h^) = 8n (2^- < 1% >) + 8n ( T (%- < T<%, >). (3.38) 
Here, the left-hand side G^vft, h^) denotes the first perturbation of the Einstein tensor 

(2) 

G^v about the background metric 7 M! , , but with its linear argument taken to be h ^ 
rather than h ( % . Thus, - 2 G (1 ^(7, h^) is given by the left-hand side of Eq.(3.20), 
subject to Eqs.(3.21,22). Hence, the left-hand side of Eq.(3.37) is linear in /r ^ and 
its derivatives, whereas the right-hand side is quadratic in first-order fluctuations. By 
contrast with Eq.(3.36), the terms in Eq.(3.37) vary over length-scales of order A . 

4. Scalar field: harmonic description 

Consider small bosonic perturbations cj)^ and /Z 1 ^ , obeying the linearised classical 
field equations (3.20) and (3.25) about a spherically-symmetric classical solution ($, 7^) 
of the Riemannian field equations (3.7-11) for Einstein gravity, coupled minimally to a 
massless scalar field. The background spherically-symmetric data for $ and 7 M „ are posed, 
as in Sees. 1,2, on the initial and final 3-dimensional boundaries, separated at spatial infinity 
by a 'Euclidean time-separation' r > . Similarly, the linearised classical perturbations 
</>W and h^liv are to be regarded as the solutions to a coupled linear elliptic problem, 
subject to prescribed linearised perturbations <f>^ (say) and h ^ on the initial and final 
boundaries. 

Because of the spherical symmetry of the background ($,7^), one may expand the 
Riemannian 4-dimensional perturbation <f>^ in the form 

1 00 l 

<P^(r,r,e,<P) = - J2 E Y ern (^)RUr,r). (4.1) 

1 = m = -£ 

Here, Ye m (fl) denotes the (£,m) scalar spherical harmonic of [48]. 

Similarly, a generic Riemannian metric perturbation ft^f may be expanded out as 
a sum over tensor, vector and scalar [l , m) harmonics, each weighted by a function of r 
and r [32,33,49-52]. The amplitudes for photon (spin-1) and graviton (spin-2) emission 
following black-hole collapse are treated in [2]; there, further details of the spin-1 and 
spin-2 harmonics are given. But note that, because of the coupled nature of the linearised 
field equations (3.20,25) for (p^ and h \ v , the resulting linear field equations for Rt m {ji r ) 
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of Eq.(4.1) and its gravitational analogues will also be coupled in the strong-field 'collapse' 
region of the (Riemannian) 'space-time'. 

The boundary conditions on the radial functions Re m (r, r) as r — > follow from the 
regularity there of the whole Riemannian solution, consisting of <p and the 4-metric g^ v 
(but viewed in 'nearly-Cartesian coordinates' near r = 0). This regularity of the solution 
in turn follows since the coupled field equations are 'elliptic modulo gauge'. For simplicity, 
the boundary data, on both the initial and the final 3-surface, should be chosen to be 
suitably regular or smooth over M 3 , in addition to being asymptotically flat. Even when 
one takes a complex Lorentzian time-separation-at-infinity 

T = rexp(-z0), (4.2) 

as in Eq.(2.3), with < 6 < tv/2 , one expects that the field equations (up to gauge) will 
be strongly elliptic [13], whence all classical fields must be analytic in the interior of the 
large cylindrical boundary formed by the initial and final surfaces, together with a surface 
at large r . 

Suppose that the boundary conditions on the final surface are taken to describe very 
weak and diffuse scalar and gravitational fields, to be regarded as perturbations of flat 3- 
space E 3 . (One also requires that the ADM mass of the final intrinsic boundary 3-metric 
gij , as computed from the 1/r part of the fall-off of gij to the flat metric dij [16,17,53], 
should be the same as the ADM mass of the initial surface. This will be discussed further 
in Sec. 5 below.) Physically, such weak and diffuse final boundary data may be imagined 
to be a possible late-time remnant of gravitational collapse, namely, a snap-shot of a large 
number of scalar particles and gravitons, as they make their way out to infinity. Near the 
final surface, the coupling in Eqs. (3.20,25) between the linearised perturbations (p^ and 
hP~\ v will almost have disappeared. The perturbed scalar field equation at late times is 
simply 

VV^ (1) = 0, (4.3) 

with respect to the spherically-symmetric background geometry 7 Mi , . 

Making the mode decomposition (4.1) of one obtains the (£ , m) mode equation 

( e ( b -^/ 2 d r ) 2 R lm + (d T ) 2 R im + l - (d T (a-b))(d T R £rn )-V i (r,r)R £rn = 0. (4.4) 
Here, the potential Vg(r, r) is given by 

W(r,r)=^(/(/ + l) + B=^), (4.5) 
and m(r, r) is defined by 

exp(-a(r,r)) = 1 - ^tlll . (4.6) 

In an exact Schwarzschild solution with no scalar field, one would have e b = e~ a = 
1 — (2M/r), with M the Schwarzschild mass; in that case, m(r, r) would be identically 
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M . The potential Vg(r, r) of Eq.(4.5) generalises the well-known massless-scalar effective 
potential in the exact Schwarzschild geometry [16], which vanishes at the event horizon 
{r = 2M} and at spatial infinity, and has a peak near {r = 3M}. 

The definition (4.6) of m(r, r) is also consistent with the usual description of the 
Lorentzian-signature Vaidya metric [19,20]. In terms of a null coordinate u and an intrinsic 
radial coordinate r , the Vaidya metric reads 

ds 2 = -2 du dr - (l - du 2 + r 2 (d9 2 + sin 2 6 # 2 ) . (4.7) 

Here, m(u) is a monotonic-decreasing smooth function of u , corresponding to a suit- 
able spherically-symmetric outflow of null particles, for example by taking the energy- 
momentum tensor of a black hole evaporating via emission of scalar particles at the speed 
of light. The Vaidya metric has been used often to give an approximate gravitational back- 
ground for black-hole evaporation at late times [54-56]. In connection with our present 
work, the Vaidya geometry has been treated in [21]. 

There is, of course, an analogous decoupled harmonic decomposition, valid near the fi- 
nal surface, for the weak gravitational-wave perturbations about the spherically-symmetric 
background — again described in [2]. For simplicity of exposition, we shall here restrict 
attention to weak-field final configurations for spin (scalar) , and calculate their quantum 
amplitudes, on the further assumption that the final 3-metric hijF is exactly spherically 
symmetric (in addition to the assumed spherical symmetry of the initial data <pi and hijj). 
Once the methods are established in the simplest spin-0 case, generalisation to the case of 
higher-spin fields becomes relatively straightforward. 

5. The classical action 

Consider, for definiteness, an asymptotically-flat Lorentzian-signature classical solu- 
tion (g^,^) of the coupled Einstein/massless-scalar field equations, between an initial 
hypersurface E/ and a final hypersurface E^ , separated by a Lorentzian proper time T 
at spatial infinity. Write 5 for the Lorentzian action functional, which corresponds to 
the Riemannian action functional / of Eq.(3.4) with suitable boundary contributions [12], 
appropriate to fixing the boundary data (h^ , 4>)i and (h^ , 4>)f , according to iS = — I . 
At the Lorentzian-signature solution above, one has [12,57] the classical action 

Sciass [(hij ,(f))r, (h^ , 4>) F ; T] 

= ^{L~D d3x * iih * + K1-/J dZx MT - (8 ' 1) 

Here, n 1 ^ = tx^ 1 is 167r times the Lorentzian momentum conjugate to the 'coordinate' 
variable hij on a space-like hypersurface, in a 3 + 1 Hamiltonian decomposition of the 
Einstein/massless-scalar theory [58]. Explicitly, in terms of the Lorentzian-signature sec- 
ond fundamental form Kij = of the hypersurface [12,41], 71^ is given by 

7r« = hS (K ij -Kh ij ), (5.2) 
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where h = det(/i^) and K = h l i K^ . Further, ir^ is the Lorentzian momentum conjugate 
to the 'coordinate' variable <j) . Explicitly, 

7T4, = hh n^V^cj) , (5.3) 

where denotes the (Lorentzian-signature) future-directed unit time-like vector normal 
to the hypersurface. 

Suppose instead that one has a complex or a Riemannian solution (g^ , <p) between 
asymptotically-flat boundary data {hij ,(/>)/ and {hij ,4>)f on initial and final hypersur- 
faces E/,Ei?, where the time-separation T at infinity has the form T = rexp(— iO), as 
in Eq.(2.3), with r positive real, (0 < 6 < n/2). As above, this is expected to provide the 
natural arena for asymptotically-flat boundary-value problems involving gravitation, pro- 
vided strong ellipticity holds, up to gauge. For such a solution, the Lorentzian-signature 
classical action S c \ a ss is still defined by Eq.(5.1). This action will in general be complex, 
although for a real Riemannian solution with 9 = tt/2 , the Riemannian action I c iass , 
defined by I c iass = — iS c \ ass , is real. The boundary contribution at spatial infinity to 
the Riemannian action functional / of Eq.(5.1) is Mr [12]. The boundary contributions 
to the functional / , due to the presence of the boundaries S j and E p with specified data 
{hij , 4>)i and {h^ , <j)) F , are 

h + I F = ^- Qf - d 3 x e w lJ h^ + ~ d 3 x e Tt^ . (5.4) 

Here, 

e7 r^ = hh ( e KV - e K (5.5) 

is given by the same formula as tc 1 ^ in Eq.(5.2), except that Kij has been replaced by the 
'Euclidean' second fundamental form e Kij , as defined and used in Eqs. (2.6.23,24) of [12]. 
In particular, 

e K i:i = -i Kij . (5.6) 

Similarly, the scalar-momentum variable tc^ of Eq.(5.3) has been replaced by its 'Euclidean' 
version e 7r^ , defined by 

e 7r^ = hh e n M V M< /> , (5.7) 

where [12] 

denotes the unit future-directed Riemannian normal. 

The quantity M in Eq.(5.1) is the ADM mass of the 'space-time', as measured near 
spatial infinity from the 1/r part of the fall-off of the intrinsic spatial metric hij on Ej 
and E_p [16,53]. As mentioned in Sec.l, it is essential, for a well-posed asymptotically-flat 
boundary-value problem, that the intrinsic metrics hijj and hijp be chosen to have the 
same value of M . Otherwise, if Mj ^ Mp , then any classical infilling 'space-time' will 
have E j and E p badly embedded near spatial infinity, and the entire 4-metric g^ v will not 
fall off to flatness at the standard 1/r rate, as r — > oo [17]. 
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(5.8) 



In applications to particle emission, including the case of nearly-spherical collapse to 
a black hole, we naturally make use of the perturbative splitting 

9»u ~ 7m- + hil iu + h {2 i v + ... , (5.9) 

rsj $ + 0« + (/,( 2 ) + ... . (5.10) 

Here, the spherically-symmetric 'background' solution (7^ , $) obeys the coupled Einstein 
/massless-scalar classical field equations, as does the full classical solution {g^ v ,</>). (The 
formal device of including a small parameter e has been relaxed here; we now set e = 1.) 

The linearised fields h ^ v and may be decomposed into sums of appropriate an- 
gular harmonics, labelled by quantum numbers (£,m), as in Sec. 4, and without loss of 
generality it may be assumed that any spherically-symmetric i = linear-order pertur- 
bation modes have been absorbed into the spherically-symmetric background (7^,$). 
Then (say) the Lorentzian classical action Sciass of Eq.(5.1) may be split as 

S 'class = Scliss + Sci^ + S^lals + • • • • 0>-H) 

Here, S^iass is the background action, given by Eq.(5.1), but evaluated for the spherically- 
symmetric solution (7^,$). The mass M appearing in S^ss wm ^ e ^ na ^ determined 

(2) 

from 'jiji or 7^ . The next term is S^iass > f° rme d quadratically from the linear-order 
perturbations; the linear-order term S^iass * s zero ? because one is perturbing around a 
classical solution. In an obvious notation, one has 




(5.12) 



Thus, S^iass is contructed only from quantities on the boundaries E/ and E_p . Note that 

(2) 

there is no contribution to the second-order expression S^ ss from the — MT term in 
Eq.(5.1), again because of the above definitions. 

The expression (5.1) for S c \ &ss [{hij ,<f>)i ; (h^ ,4>)f', T], together with the asymptotic 

(2) 

series (5.9) for the classical action and the expression (5.10) for S^ ss formed from the 
linearised perturbations, will be basic in calculations concerning quantum amplitudes in 
subsequent work. 

6. Adiabatic radial functions 

We return to the evolution of linearised scalar-field perturbations <jft-\ following the 
mode sum (angular decomposition) of Eq.(4.1). For the quantum amplitudes of interest, 
we must compute expressions of the form 

Amplitude = const, x expjz [{hij ,4>) I ; [h^ ,4>) F ; T]|, (6.1) 

where, equivalently, i<S c iass = — /class • As above, the time-interval T , measured at spatial 
infinity, must be of the form T = \T\ exp(— 16), (0 < 9 < n/2), for a classical solution 
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to exist for the boundary-value problem. The classical action S^iass is given in Eq.(5.1) 
in terms of integrals taken over the boundaries £/ and T, F , subject to the classical field 
equations. 

Consider the amplitude corresponding to weak-field non-spherical data {hP^)F and 

(4>^)f ° n the final surface, given at lowest order by exp(i5^ ss ) . For simplicity, take the 
initial data to be exactly spherically symmetric, namely, (7^ , $)/ . Equivalently, 

h { % = ; </> (1 ) = 0. (6.2) 

The amplitude exp(i5' c i ass ) will then depend only on the contributions at the final surface 
Hp in Eq.(5.12) [which themselves depend on {h \j F , 0^ ; T)]. As a practical matter, 

one could easily put non-zero (h , </> ( - 1 - ) )/ back into the calculations that follow. Phys- 
ically, the analogous step of 'turning back on the early-time perturbations' corresponds, 
in 'particle language' rather than in the 'field language' being used in this paper, to the 
inclusion of extra particles in the in-states, together with the original spherical collapsing 
matter, and asking for the late-time consequences. This was first carried out by Wald [59]. 

In this paper, we concentrate only on the scalar-field contribution to the quantum 
amplitude exp(i5' c i ass ). That is, we compute 

Sil„ca,„ = U <?* *?> (6-3) 

where the linearised perturbations (/Z 1 ^ ,4>^) obey the linearised field equations (3.20- 
22,25) about the spherically-symmetric background (7^,$). Here, (h^v , <j>^) must 

agree with the prescribed final data {h^l- at the final surface £f , and be zero 

at the initial surface £/ . In the Riemannian case, with a real Euclidean time-interval 
r between £/ and £f , or in the case (2.3) of a complex time-interval T = rexp(— z#) 
between £/ and £i? , with < 9 < tt/2 , this linear boundary-value problem is expected 
to be well-posed. The other, gravitational, contribution 

S {2) = — / d 3 x tt (1) ^' h {1) (6 4) 

to S^ ss in Eq.(5.12) is studied in [2]. 

Following Sec. 4, at late times the perturbed scalar-field equation reduces to 

V^V^ (1) = 0, (6.5) 

with respect to the spherically-symmetric background 7 Mi , . Here, in contrast to Eq.(3.5), 
it is more suitable to work with the Lorentzian background gravitational field: 

ds 2 = _ e Ht,r) dt 2 + e a(t,r) dj ,2 + r 2 ^2 + ^2 Q (g g) 
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By analogy with Eq.(4.1), one makes the mode decomposition with respect to 'Lorentzian 
coordinates' {t,r,6,<fi): 

oo £ 

<P^(t,r,e,<P) = - E Y em (Q)R em (t,r). (6.7) 

1 = m = -£ 

As in Eq.(4.4), one arrives at the (£ ,m) mode equation: 

{e {h - a ^ 2 d r ) 2 R, m - (d t ) 2 R £m - l(dt(a-b))(d t R tm )-V e (t,r)R tm = 0. (6.8) 

Here, Vi(t,r) is defined by Eq.(4.5), except that, in all its appearances, the argument r is 
replaced by t . Similarly, the function m(t,r) is defined by Eq.(4.6), with a corresponding 
replacement of r by t . 

For high frequencies of oscillation in the nearly-Lorentzian case, with small angle 9 
of rotation into the complex, it becomes simpler to understand the solutions of the mode 
equation (6.8). Consider a solution Re m (t, r) of the form 

Remit, r) ~ exp(ifct) f fc £ m (t,r), (6.9) 

where k is a 'large' frequency, but where, in contrast, £k£m(t, r) varies 'slowly' with respect 
to t . In particular, we require that, near spatial infinity, with r — > oo , R£ m (t, r) reduces 
to a flat space-time separated solution, in which £kem(t,r) loses its t-dependence [see 
Eqs.(6.11,14) below]. 

Our boundary-value problem is for scalar perturbations <j>^(t, r, 9, 0), or equivalently 
for functions Ri m {t,r) as in Eqs.(6.7,8), subject to the initial condition (j)^'\ t _ Q = and 

to prescribed real final data (f>^\ t=T ■ Were the propagation simply in flat space-time, 
the solution would be of the form 

= - E E r dk a ^m tkem(r) Y im (Sl), (6.10) 

where the {ak£ m } are real coefficients and each function £,k£m( r ) is proportional (up to a 
factor of r) to a spherical Bessel function je(kr) [60]. In our gravitational-collapse case, 
ik£m becomes a function of t as well as of r , but otherwise the pattern remains: 

= - E E f°° dk ak ^ ^£m(t,r) Y tm (Cl). (6.11) 

r £ = m = -£ J ~°° Sm{Hl) 

Here, the {cikem} characterise the final data: they can be constructed from the given 
4>^\ t _ T by inverting Eq.(6.11). The functions £,k£m(t,r) are defined in the adiabatic or 
large— | A;j limit, as in the previous paragraph, via Eq.(6.9), where Rg m (t,r) obeys the mode 
equation (6.8). 
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More precisely, provided that k is large, in the sense that the adiabatic approximation 



\k\ > - |d - b\, 



(6.12) 



\k\ > 



skim 



klm 



k 2 > 



Ikim 



holds, the mode equation reduces approximately to 



or V or / 



(6.13) 



(6.14) 



Of course, the functions e^ b a ^ 2 and Ve do still vary with the time-coordinate t , but only 
adiabatically or 'slowly'. 

As described further in [1,5,19], the geometry in the radiative region of the space-time 
is expected to be approximated very accurately by a spherically-symmetric Vaidya metric 
[19,20], corresponding to a luminosity in the radiated particles which varies slowly with 
time. Such a metric can be put in the diagonal form 



e~ a = 1 



2m(t, r) 



e b = 



m 



(6.15) 



\k\ > 


m 






m 



where m(t,r) is a slowly-varying function, with m = (dm/dt), and where the function 
f(m) depends on the details of the radiation. Then Eq.(6.12) implies that 



(6.16) 



provided that 2m(t, r) < r < 4m(t,r). In this case, the rate of change of the metric 
with time is slow compared to the typical frequencies of the radiation; further, the time- 
variation scale of the background space-time metric 7 Mi , is much greater than the period 
of the waves. With frequencies of magnitudes \k\ ~ m _1 dominating the radiation, and 
with 1 771 1 of order m~ 2 [61], the adiabatic approximation is equivalent to m 2 ^> 1 , which 
corresponds to the semi-classical approximation. If, as expected [61], m 3 is a measure of 
the time taken by the hole to evaporate, then r < 4m m 3 , provided that m 2 ^> 1 . 
Thus, in the large-/c approximation used in deriving Eq.(6.14), it is valid at lowest order 
to neglect time-derivatives of the background metric, out to radii small compared with the 
evaporation time of the hole and with the time since the hole was formed. 

It is natural to define a generalisation r* of the standard Regge-Wheeler coordinate 



r* s for the Schwarzschild geometry [16,49], according to 



d 

dr* 



= e (b-a)/2 

dr 



(6.17) 



Under the above conditions, the time-dependence of r*(t,r) is negligibly small, and one 

has r* ~ r* 



s for large r , where, by definition, 

r* s = r + 2Mlog((r/2M) - l) 



(6.18) 
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is the Regge- Wheeler coordinate, expressed in terms of the Schwarzschild radial coordinate 
r. In terms of the variable r* , the approximate (adiabatic) mode equation (6.14) reads 



%^ + (k 2 -V i )Z kem = 0. (6.19) 

7. Boundary conditions 

We now consider, in more detail, a set of suitable radial functions {iktmir)} on the 
final surface Tip . As above, since the mode equation (6.8) does not depend on the quantum 
number m , we may choose $,k£m( r ) = £,ke( r ), independently of m . 

We seek a complete set, such that any smooth perturbation field <p^\T, r, 9, <j>) of 
rapid decay near spatial infinity, when restricted to the final surface {t = T}, can be 
expanded in terms of the iktmir). The 'left' boundary condition on the radial functions 
{£,ke(r)} is that of regularity at the origin r = : 

6w(0) = 0. (7.1) 

The solution to the radial equation, regular near the origin, is: 

Zkt(r) = r<p ke (r) oc r j e (kr) oc { (const, x {kr) t+1 ) + 0((kr) e+3 ) }, (7.2) 

where, again, jg denotes a spherical Bessel function [60]; we have assumed that, for small 
r , one has m(r) oc r 3 , and we have neglected 0(r 2 ) terms. These radial functions are 
purely real, for real k and r . For k purely real and positive, the radial functions describe 
standing waves, which, for mode time-dependence e ±lkt , have equal amounts of 'ingoing' 
and 'outgoing' radiation. 

For the 'right' boundary condition, note that the potential Vi(r), following from 
Eq.(4.5), vanishes sufficiently rapidly as r — > oo that a real solution to Eq.(6.19) obeys 

iktir) ~ (z ke exp(z£r* s ) + z* M exp(- z/cr* s ) j (7.3) 

as r — > oo . Here, the z k e are certain dimensionless complex coefficients, which can be 
determined via the differential equation by using the regularity at r = . The (approxi- 
mately) conserved Wronskian for Eq.(6.19), together with Eq.(7.3), and the property 

exp(z(A; — k') r* ) , 
Jim P 7, ' s) = iv5{k-k% (7.4) 

give the normalisation condition, for — oo < k , k' < oo and Rqo — > oo : 

dr e^' 2 i u [v) H t {T) = 2tt \z m \ 2 (6(k - k') + S(k + k')). (7.5) 



I 

Jo 



19 



This normalisation is only possible within the adiabatic approximation. Note that the 
radial functions {ikt} form a complete set only for k > , as a result of our boundary 
conditions. 



bution to the total classical Lorentzian action S c \ ass = S^iass + ^cfass + • • • °^ Eq.(5.11), 

(2) (2) 

with S^iass gi ven by Eq.(5.12). This contribution, namely S^iass scalar of Eq.(6.3), is given 



The above result makes it possible to evaluate the perturbative massless-scalar contri- 
on to the total classical Lc 

iS™. given by Eq.(5.12). 
in the notation of Eq.(6.7) by 

; T] = \ E E / ^ e(a_b)/2 I T > (™) 



since 



y ^ Y £m Y[, m , = Sw $mm' ■ (7.7) 

Within the adiabatic approximation above, and using Eq.(7.5), this gives the frequency- 
space form of the classical action: 

s i\Ls[{ a kem} ] T] = 7T E E / dk k\z k e\ \a k e m + a- k £m\ cot(fcT), (7.8) 



in terms of the final data {a^ m }. 

From a mathematical point of view, one would expect to work only with the set 
of square-integrable scalar wave-functions on the final boundary Hp , that is, the set 
L 2 (R 3 , dref"-^ 2 ). To express this, define 

^m(r) = rjdQ Y lm {Q) ^\t,r^)\ t=T . (7.9) 



Then the square-integrability condition reads 

V / 

2tt 



£m 

or, equivalent ly, 



1 fRoo 

r E / dr e^/ 2 \^ m (r)\ < oo , (7.10) 



/oo 
dfc |z^| |a fc £ m + a_ fc £ m | < oo . (7.11) 



The left-hand sides of Eqs.(7.10,ll) are in fact equal. This arises from the complete- 
ness property 



( a-6)/2 r dk Hkt(r) H k t(r') = 4n6{r _ r/) (712) 

J-oo \Zk(\ 
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and the inverse of Eq.(6.11): 



^ /'•too 

a kim + a- kim = ■ / * e (a " 6)/2 ^(r)^ m (r). (7.13) 

27v\z k e\ Jo 

From a physical point of view, one expects also that taking scalar boundary data 
which are not square-integrable will lead to various undesirable properties, such as infinite 
total energy of the system, or an infinite or ill-defined action. 

8. Analytic continuation 

The perturbative classical scalar action S^iass of Eq.(7.8) was derived subject to the 
adiabatic approximation and to the requirement that the time-interval T between the 
initial and final surfaces, measured at spatial infinity, should be complex, of the form 
T = |T|exp(— i9), with < 9 < n/2 . In this case, the term kcot(kT) in the integrand 
of Eq.(7.8) remains bounded near k = 0, and one expects to obtain a finite complex- 
valued action S^l ss [{akem} ', T], given square-integrable data on the final surface • 

(2) 

Further, the dependence of the complex function S^ aa [{akem} ] T] on the complex variable 
T is expected to be complex-analytic in this domain (0 < 9 < 7r/2), and, following 
Feynman [14,15], ordinary Lorentzian-signature quantum amplitudes should be given by 

( 2) 

the limiting behaviour of exp^S^^J as 9 — > 0+ . 

If, on the other hand, one restricts attention to the exactly Lorentzian-signature case 
(9 = 0), then the integral in Eq.(7.8) will typically diverge, due to the simple poles on the 
real-frequency axis at 

(n = 1,2,3, ... ). (8.1) 

We now restrict attention to the case in which T is only slightly complex, writing 
T = \T\ exp(— i5), with < 5 1 . The spherically-symmetric 'background' 4-geometry 
7^ and scalar field $ will then be complex. Consider an integral such as Eq.(7.8) for 

4ass[K^}; T ]- Write this as 

poo 

J = / dk f em (k)cot(kT), (8.2) 

im J ° 

where 

I 1 2 i 1 2 

fim(k) = 7ik\z k e\ \akim + a-Um\ ■ (8.3) 

There are infinitely many simple poles of the integrand at k = k n (n = 1, 2, . . . ), just 
above the positive real fc-axis. We then deform the original contour C along the positive 
real /c-axis into three parts, C e , Cr and C a , where < a < 1. The contour C e lies in 
the lower half-plane, half-encircling each of the simple poles near the positive real fc-axis, 



k — k n — — 
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with radius e . The curve Cr , also in the lower half-plane, is an arc of a circle \k\ = R of 
large radius. The curve C a is part of the radial line arg(/c) = — a . We write 



'c a +c R -c e (8.4) 



^ Jc a +c R -c r 

= J a + Jr + J e . 

Starting with the integral Jr , one finds 



\Jr\ < J2 / M R\fem(R,d)\ coth(\T\R 8m6), (8.5) 
em Jo 

where k = Re~ lB on Cr , and we have used |cot(fcT)| < coth(|T|i?sin#). One expects 
that, when the limit R — > oo is eventually taken, the contribution from Cr to the total 
action should vanish; this requires that \fim{k)\ should decay at least as rapidly as \k\ , 
as \k\ — )■ oo . In fact, on dimensional grounds, one expects that 

\fim{k)\ ~ const, x |/c|~ 3 (8.6) 
as \k\ — > oo . To see this, rewrite the radial equation (6.19) in terms of the operator 

£ < = e{h ~ a, ' 2 l{ e """ ,l2 i^)- v ' (r) ' (8 - 7) 

which is self-adjoint with respect to the inner product in Eq.(7.5). Then note that Eq.(7.13) 
can be rewritten as 

a k e m + a- k e m = ■ ^ / dr e (a_&)/2 £ ke {r) C e ip im (r). (8.8) 

2nk 2 \ zu\ Jo 

We have used the boundary condition (7.1) and assumed that ip£ m (r) dies out at large 
r . The form (8.8) is just an expression of the self-adjointness of the radial equation. 
Now consider the dimensions of the quantities involved. One finds [1] that ip£ m (r) has 
dimensions of length and that \Zkt\ is dimensionless. In the limit Roo — > oo , and for 
large k (so taking a WKB approximation for the radial functions), the integral in Eq.(8.8) 
can only involve the dimensionless frequency 2Mk , where M is the total mass (true ADM 
mass) of the space-time. This gives the desired behaviour (8.6) at large \k\ . 

The contour C e gives a purely imaginary contribution to the total Lorentzian action; 
also (see below), the curve C a gives a complex contribution. We shall interpret the quan- 
tity exp[— 2 Im(S')], up to normalisation, as describing the conditional probability density 
over the final boundary data. To compute J e , we assume that f(. m {k) is analytic in a 
neighbourhood of k = a n , where 

Tl 7T 

a n = — (n = 1,2,3, ...). (8.9) 
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Note the difference between the definitions (8.1) of k n and (8.9) of a n . Then 

J r = - lim V / dk f £m (k) cot(k\T\) 

6 > ^ I r~i 

0™ ^ O e 



\T 

For the curve C a , one has 



e 

f ----- • 

em 

oo 

= W\ E E ft™(Vn)- 
£m n = l 



(8.10) 



J Q = - V / e- lQ /, m (|A;|,a) cot(|fc| e" lQ |T|) . (8.11) 
We shall need the properties [60] 

V 

n = — oo v ' 

and 

= P.P.- =f iiv5(x-a), (8.13) 



(x — a ± ze) (x — a) 

where P.P. denotes the principal part. Assuming that ft m (k) is regular along C a , one 
has, for small a : 

■#>= (l-ia)£ E /Vl J:f kla \ x (8-14) 

a ^0+ Jq /cT — 77-7T — ZCK 

n = — oo w V 1 / 

In the further limit R — > oo , this gives 

oo 

J a = P.V. + ^ E (8-15) 



£m n = 1 



where P.V. denotes the principal-value part of the integral. 

Using Eqs.(8.4,5,10,15), the classical action for massless scalar-field perturbations, 
when T = \T\ exp(— id) is very slightly complex, is 

2 • oo t oo 

S %L[i a k£m}; \ T \] = real P art + -j^r E E E 

I I ^ = m=-i n = l 



(8.16) 

= real part + — — ^ °n |aw m + a_ n £ m | . 



2*7T 2 



Imn 



The real part of S^fass is, of course, also calculable from the equations above. It con- 
tains the principal-value term and the real part of Eq.(8.10). The main, semi-classical, 
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contribution to the quantum amplitude is then exp(iS£ ls ^ [{afe£ m } ; \T\]). The probability 

(2) 

distribution for final configurations involves only Im(S'^ lass ); the more probable configu- 

(2) 

rations have S^ ss lying only infinitesimally in the upper half-plane. Whether probable 
or not, those final configurations {akim} which contribute to the probability distribution 
must yield finite expressions in the infinite sums over nl in Eq.(8.16). There will be a 
corresponding restriction when the data are instead described in terms of the spatial con- 
figurations {i()£m(r)}. Also, as can be seen in [5], the complex quantities z n i (aw m + a_ n ^ m ) 
appearing in Eq.(8.16) are related to Bogoliubov transformations between initial and final 
states, thus providing a further characterisation of the finiteness of Im^j^J in Eq.(8.16). 

With regard to the sum over i in Eq.(8.16), one imagines that a cut-off ^ max can be 
provided by the radial equation (6.19). In the region where (Vg(r) — k 2 ) > 0, one has 
exponentially growing radial functions, whereas for (Vg(r) — k 2 ) < one has oscillatory 
radial functions. One defines £ max by (V£ ma , x (r) — k 2 ) = and restricts attention mainly 
to oscillatory solutions. 

When one has both initial and final non-zero Dirichlet data labelled by 'coordinates' 
{ a< ktm} an d { a ifoL}' the perturbative classical scalar action S^ ss includes separate terms 
of the form (8.16) for the initial and final data. But S^l ss also includes a cross-term 
between a^ m and , which represents the correlation or mixing between the initial 

and final data. The total action will naturally be symmetric in and o/y^ , and the 
coefficients z n g will be the same (they are time-independent) up to a phase. For large 
\T\ , the cross-term becomes negligible, and one has two independent contributions to the 
classical action, one being a functional of {o^ m }, the other of {a^ml- 

9. Conclusion 

In this paper, we have derived through Eq.(8.16) the quantum amplitude for a spheric- 
ally-symmetric configuration (hij ,(/>)/ on the initial surface E/ to become a configuration 
(hij,(p)F on the final surface T,p , with Lorentzian time-interval T at spatial infinity. 
Here, 4>p will, in general, be anisotropic, although (for simplicity) we assumed that the 
final 3-dimensional metric hijp is also spherically symmetric. In the amplitude, which 

for a locally-supersymmetric theory is proportional to exp^S'^J, the classical action 
depends approximately quadratically on the (non-spherical) perturbative part of the final 

(2) 

data 4>f ■ Further, S^ ss has both a real and an imaginary part. The imaginary part leads 

2 (2) 

to a Gaussian probability density |$| oc exp(— 2 Im^^J), while the real part gives rapid 
oscillations in the phase of the quantum amplitude or wave function $ . 

We have arrived at a quantum amplitude (not just a probability distribution) for 
such processes, simply by following Feynman's +ie prescription, applied to the exactly 
semi-classical expression (2.2) for the quantum amplitude. This, in turn, is derived via 
Dirac's canonical-quantisation approach, for a locally-supersymmetric Lagrangian such 
as that of gauge-invariant N = 1 supergravity. The boundary conditions are treated by 
rotating the time-interval T into the lower-half complex plane: T — > |T| exp(— iO), for 
< 6 < n / 2 . We then studied the classical and corresponding quantum-mechanical 
boundary-value problems, before rotating 6 back towards zero. 
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These ideas have also been applied to black-hole evaporation for particles of spin 1 and 
2 [2], and to the fermionic spin-| case [3]. But, in those references, the form of the complex 
quantum amplitudes was not derived in the greater detail given in the full computation of 
spin-0 amplitudes of the present paper. In [5,6], we made a connection relating the present 
description and calculation of quantum amplitudes to the familiar description in terms of 
Bogoliubov coefficients [62-64]. A more general conceptual framework has been provided 
within the language of coherent and squeezed states [65,66]. 
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